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ABSTRACT 



Aims. By measuring the gravitational redshift of spectral features in the spectrum of thermal radiation emitted from neutron stars, use- 
ful constraints for the equation of state of superdense matter can be obtained via an estimate of the mass-to-radius ratio. We searched 
for spectral line features in the high-resolution X-ray spectrum of the isolated neutron star RX J0720.4-3125. 
Methods. Our target was observed by XMM-Newton on many occasions. We used the XMM-Newton SAS task rgscombine in order to 
create a co-added RGS spectrum of RX J0720.4-3125. We modeled the resulting spectrum with absorbed blackbody radiation with 
a Gaussian absorption features using the XSPEC package. 

Results. We found a narrow absorption feature at 0.57 keV in the co-added RGS spectrum of the isolated neutron star RX 
J0720.4-3125 with an equivalent width of 1.35 ± 0.3 eV and FWHM ~ 6.0 eV. The feature was identified with an absorption 
line of highly ionized oxygen O VII, most probably originating in the ambient medium of RX J0720.4-3125. An extensive investi- 
gation with the photo-ionization code CLOUDY indicates the possibility that the optical flux excess observed in the spectrum of RX 
J0720.4-3125 at least partially originates in a relatively dense (e.g. iiu ~ lO^cm"^) slab, located in the vicinity of the neutron star 
(e.g. ~ lO'^cm). 
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1. Introduction 

The study of thermally emitting and radio-quiet, nearby, isolated 
■ neutron stars (INSs) may allow an important input to our under- 
standing of neutron stars. In particular, measuring the gravita- 
tional redshift of a spectral feature in the spectrum of thermal ra- 
diation emitted from the neutron star surface (atmosphere) may 
provide a useful constraint for various equations of state for su- 
perdense m atter. Independent of the estimate of the neutron star 
radius (e.g. Triimper et al.l 12004) from the thermal spectrum of 
an INS, it will also allow us to directly estimate the mass-to- 
' radius ratio. RX J0720.4-3125, originally discov ered as a soft 
X-ray source during the ROSAT All-Sky Survey bvl Haberl et alj 
( Il997h . is a special case among the INS. Clear spectral vari- 
ations on times scales of years were detected during XMM- 
Newton observations in the data of the high s pectral resolution 
Reflection Grating Spectrometers (RGS) (d e Vries et al.ll2004 



Table 1. XMM-Newton 
J0720.4-3125 



RGS 1/2 observations of RX 



and th e imaging-spectroscopic instrument EPIC pn (i Haberl et alj 
|2006|) . However, during individual RGS observations the num- 
ber of counts received from RX J0720.4-3125 in a single reso- 
lution element were not enough to study the spectrum in detail. 

From this point of view, it is interesting to study a co-added 
XMM-Newton RGS spectrum of RX J0720.4-3125, despite of 
the spectral variations shown by the INS. 



Obs. ID 


Obs. Date Start 


Exposure 


Effective exposure 




JD - 2400000 


ksec 


ksec 


0124100101 


51678.1503935 


65.87 


15.45 


0132520301 


51870.2714931 


30.91 


30.24 


0156960201 


52585.2443171 


30.24 


29.60 


0156960401 


52587.3090509 


32.04 


31.03 


0158360201 


52762.0499190 


81.63 


73.43 


0161960201 


52940.2405093 


44.92 


44.67 


0164560501 


53147.9273380 


51.95 


44.36 


0300520201 


53488.8618634 


53.31 


50.69 


0300520301 


53636.4891551 


53.01 


51.38 


0311590101 


53687.4349306 


39.71 


38.91 


0400140301 


53877.6977662 


21.91 


21.57 


0400140401 


54044.9718634 


21.91 


21.69 


0502710201 


54226.2093171 


21.91 


21.57 


0502710301 


54421.7184259 


24.92 


24.59 



* Based on observations obtained with XMM-Newton, an ESA sci- 
ence mission with instruments and contributions directly funded by 
ESA Member States and NASA 



2. Observations and data reduction 

RX J07 20.4-3125 has be en observed many times by XMM- 
Newton toohleet a l. 2009) and we focus here on the data col- 
lected with RGS fden Herder etal.ll200lh from the 14 publicly 
available observations, in total presenting about 500 ks of effec- 
tive exposure time. 

The data were reduced using standard threads from the 
XMM-Newton data analysis package SAS version 8.0.1. We re- 
processed all publicly available data (see Table [T]) with the stan- 
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Fig. 1. XMM-Newton high-resolution co-added spectrum of the 
isolated neutron star RX J0720.4-3 125. An absorbed blackbody 
(blue) and an absorbed blackbody with two additional Gaussian 
absorption lines (red, a wide absorption feature at 0.3 keV and 
a narrow feature at 0.57 keV) are shown. The equivalent width 
and the FWHM of the narrow absorption line (green), which we 
associate with OVII, were determined at 1.35 eV and ~6.0 eV, 
respectively. Different spectral line positions at gravitational red- 
shift (gr - 0) are indicated. 



dard metatask rgsproc with the aim to co-add the spectra with 
the task rgscombine. To determine good time intervals free of 
background flares, we used the filtering expression 'rate < 1 .0' 
on the background light curves. This reduced the total exposure 
time by ~ 13%. First order source and background spectra were 
produced from the cleaned events. A model background spec- 
trum was also generated with the task rgsbkgmodel that also was 
used for the background correction as an additional check for 
possible phenomena related to background variations seen dur- 
ing some of the observations. No significant differences between 
RGS spectra reduced in these two different ways were seen. 

The background subtracted co-added spectra of RX 
J0720.4-3125 revealed the existence of a relatively narrow 
absorption-like feature at ~0.57 keV. No sharp variation is 
present in the spectral response of the instrument, indicating an 
astronomical origin. A similar narrow absorption feature was re- 
ported in the XMM-Newton R GS spectrum of ano t her IN S RX 
J1605.3-H3249 at 0.576 keV bv lvan Kerkwiik et al.1 (l2004l) . 

We fitted the resulting co-added RGS 10 spectrum of RX 
J0720.4-3125 (see Fig. [B using the XSPEC package ver- 
sion 12.5.0. We fitted an absorbed blackbody model with a broad 
(0.08 keV) Gaussian absorption feature at 0. 3 keV, \yell kn own 
fxomXMM-Newton EPIC pn spectral studies ( lHaberlll2007l) . 

Another additive Gaussian absorption component was used 
to fit the feature at 0.57 keV. We used this model to estimate the 
equivalent width of the detected absorption feature at 0.57 keV. 
As the continuum level the absorbed blackbody model including 
the broad Gaussian line in absorption was used. With different 
levels (normalizations) of continuum and absorption lines an er- 
ror of the equivalent width was determined. Results of the fit and 
equivalent width estimates are presented in Fig.[T] 



' The RGS2 spectrum was used only for consistency checks; in the 
energy range of 0.5-0.59 keV it has a significantly reduced effective area 
owing to the failure of CCD4 in September 2000, see XMM-Newton 
Users Handbook 
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Fig. 2. Equivalent width of the absorption feature at 0.57 keV 
and blackbody temperature derived from the RGS spectra of RX 
J0720.4-3125 from three different epochs. 

As already mentioned, RX J0720.4-3125 showed spectral 
variations, first noticed from the comparison of RGS spectra ob- 
served at diff'erent epochs (d e Vries et al.ll2004l) . After the detec- 
tion of the absorption feature at 0.57 keV in the co-added RGS 
spectrum of RX J0720.4-3125, we examined its behavior over 
time. We divided all publicly available observations into three 
groups from different epochs. In the first group, we combined the 
RGS spectra of RX J0720.4-3125 from the year 2000 (Obslds 
0124100101 and 0132520301), for the second group the spectra 
between Nov. 2002 and May 2003 (0156960201, 0156960401 
and 0158360201, and in the third group, the spectra from Oct. 
2003 to Nov. 2007. This devision was inferred from our re- 
cent study on t he nature of g lobal variations observed from RX 
J0720.4-3125 dHohle et al .1 2009) which was mainly based on 
the analysis of the EPIC pn data. In this way the variations in 
the blackbody temperature and the depth of the broad absorp- 
tion feature in the spectra wi thin each group of observations 
(see. Fig. 7, ' Hohle etaTI 12009 *) are kept to a minimum and be- 
low the sensitivity of the RGS instrument. The results of this 
analysis are presented in Fig|2] It shows a tendency for a recip- 
rocal dependence of the measured equivalent width on the de- 
rived blackbody temperature. To decide whether this is due to 
pure flux (blackbody temperature) variations or to changes in 
the number of atoms in the line of sight, we fitted the three data 
sets with linked parameters of absorption features and the black- 
body temperature was left free. The estimated equivalent widths 
are consistent with the variations when parameters of the narrow 
absorption feature were left free (Fig. |2]i. However, this needs 
confirmation with higher S/N ratio data sets. 

3. Discussion 

We report the detection of an absorption line at 0.57 keV in the 
XMM-Newton RGS spectra of RX J0720.4-3125. Since there 
is no uncalibrated instrumental feature in the energy range of 
0.55-059 keV (21-22. 5A), this line is clearly of astronomical 
origin. Whether it originates at the surface (atmosphere) of the 
neutron star, in its vicinity or in the interstellar medium is not 
clear a priori. 

According to the atomic database ATOMD^ a few rela- 
tively strong spectral lines are expected in the energy range of 
0.56- 1 .0 ke V (with a typical value of gravitional redshift of INS 

" http://cxc.harvard.edu/atomdb 
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~ 0.2) from plasma with a temperature of ~ 10^ K, mainly 
belonging to highly ionized oxygen (OVII, OVIII) and neon 
(NelX). Among them the OVII resonance K,,, spectral line at 
0.574 keV has a maximal intensity and is compatible with our 
detection without any significant redshift. It is likely that the 
absorption feature consists of a number of blended absorption 
lines. However, given the data quality (i.e., spectral resolution 
and S/N ratio) it is difficult to identify all of them. 

Moreover, similar absorption features were detected in high- 
resolution X-ray spectra of other galactic a nd extragalacti c 
sources (but not in a l l case^ see, eg., Futa moto et alJ ( 2004 ). 
Nicastro et al. (2005), Yao & Wang! 120051) . iFang et alj hood) . 
Bregman & Llovd-Davies (2007B). This, as well as its rela- 
tively narrow nature dChang et al.ll200d : lPaerels|[T997h . make it 
more plausible that it originates in the circumstellar/interstellar 
medium along the line of sight towards RX J0720.4-3125. 
Indeed, the estimates of strength of the magnetic field of B,/p ~ 
2.4 X 10'^ G (from dipole braking) and ~ 5.6x'^ G (proton- 
cyclotron absorption feature) of RX J0720.4-3125 predict 20- 
60eV broadening of spectral features originating in the magne- 
tized atmosphere of the INS (Paerels 1997). Nevertheless, with 
a gravitational redshift of INS of ~ 0.16, a strong lines of OVIII 
at 0.653 keV are expected at 0.56 keV. In this case the OVII 
resonance K(j, spectral line from the INS atmosphere is ex- 
pected at 0.49 keV. A very faint absorption feature is present in 
the spectrum of RX J0720.4-3125(Fig.[B. However, given the 
fitted blackbody temperature (~ IMK) and the magnetic field 
strength (significantly increasing the binding energy of atoms) 
of RX J0720.4-3125, it is quite difficult to reconcile the pres- 
ence of significant number of O VIII or other ions (e.g., Ne IX, 
with higher excitation temperature) in the atmosphere of RX 
J0720.4-3125. 

If the detected absorption line has an origin in the cir- 
cumstellar/interstellar medium, it may be identified with the 
highly ionized oxygen O VII. Indeed, according to the atomic 
database (ATOMDB) the strongest laboratory-measured spectral 
lines from O VII in the energy range 0.55-0.59 keV are: 0.561 1, 
0.5686 and 0.5740 keV. However, the first two spectral lines are 
forbidden and only in specific conditions (density ^10'" cm"^, 
as well as temperature ~ 1.5 MK of "hybrid plasma", i.e. photo- 
ionized and collisionally excit ed) may hav e a comparable inten- 
sity to the K^resonance line dPorquet & D ubau 2000, see their 
Fig. 11). Therefore, these two lines also may be excluded from 
further consideration, unless there is evidence of the presence of 
nearby, very high-density plasma of RX J0720.4-3 125. 

From the measured value of the equivalent width We - 
1.35 + 0.3 eV of the spectral line of OVII we estimate the 
number of ions (column density, Nqvii) in the direction of RX 
J0720.4-3 125 in two different ways. First, using the dependence 
of the equivalent width of the absorptio n feature on the column 
density of OVII (Futamot o et al J 120041) . we estimated a lower 
limit of log No VII = 15.8 - 16.2 cm~^ in the direction of RX 
J0720.4-3125. 

Second, we arrive at similar results, log Nqvii ~ 15.9 - 
16.1cm"^, using the general relationship between measured 
equivalent width and number of atoms in the line of sight, 
WxIA = 8.85 X NjAfjk, where fjk = 0.695 is the oscillator 
strength of the OVII K„hne (ATOMDB) and N^ is the column 
density (ISpitzeij|l978l) . 




Fig. 3. Ionization parameter ^ - LxKnnr ) of oxygen for differ- 
ent nH= 10' - lO'cm"-' volume densities located at the distance 



of r = 10"*cm from ionizing luminosity Lx 
thermally emitting (T ~ lO^K) INS. 



10 erg s of a 



On the other hand, from the the averagcl novii= (1.35 - 
2.84) X 10 ;^cm-3 and the dis tance to RX J0720.4-3125 (d = 
360+1™ pc. lKaplan et"ani2007h one obtains a value of 15.0-15.7 
cm"^ for log Nqvii- Moreover, using a very recent estimate of 
the spatial density of the hot gas on t he Galactic plane ( nH= 
(0.3 - 3.4) X 10~^cm-^ jYao et al.l (l2009l )) based on 5uzaku obser- 
vations, we have the value of 14.0-15.3 cm"^ for log Noxygen in 
the direction of RX J0720.4-3125. 

This difference may be due to the origin of the absorption 
lines in the immediate, ambient environment of the INS. If the 
absorption feature at 0.57 keV originates in the ionized plasma 
near the INS and belongs to highly ionized O VII, then the pres- 
ence of the absorption line itself may constrain the possible loca- 
tion (separation fro m RX J0720 . 4-312 5) of the ionized cloud. 
Indeed, according to lTarter et al.l (1 19691) . if the plasma is located 
at a distance of r from RX J0720.4-3125, the ionization param- 
eter ^ - Lxl(nHr^) - LxKNnr) is strongly determined by the 
luminosity of ~ lO^^ergs"'. If s o, oxygen will be fully ph oto- 
ionized at a distance r < 10 cm (iKalknan & McCravllT982l see, 
also Fig.[3]l. 

Moreover, in order to have a sufficient amount of O VII ions 
near a thermally emitting ionizing source with temperature ~ 
lO^'K and X-ray luminosity of Lx ~ lO-'^ - lO-'^'erg s"', for the 
ionizing parameter the condition 1 < ^ < 10^* should be satisfied 
(Fig. [3j. This means that a cloud/slab with rather high density 
(eg. nn- 10' - lO^cm"-') must be located within a distance of 
10'" < r < lO'^cm from the ionizing source. 

On the other hand, the hydrogen column density must be 
< 2 X 10^" cm"^ (from the EPIC pn spectroscopy), ie. assum- 
ing oxygen abundance relative to the hydrogen n(0)/n(H) = 
4 X 10"** (Anders & Grevesse 1989), we obtain for oxygen a 
total column density A^oxygen < 8 x 10'^ cm"^. 

Given the abovementioned general restrictions, a photo- 
ionized plasma associated with RX J0720.4-3125 may con- 
tribute to the observed fluxes in the UV and optical bands. 

The optical and ultraviolet fluxes observed from rather 
faint, optical counterparts of thermally emitting INSs, in gen- 
eral, lie a factor of ~ 1 above the e xtrapolated blackbody 
spectrum in X-rays (e.g. iKaplanI l2008h . the so-called "opti- 
cal excess". In the case of RX J1856.4-3754, the optical/UV 



3 IFutamoto et al.l | |2004|) did not detect th e absorption line of O VII 
in the spectrum of Cyg X-2. According to lYao & Wan^ d2005h . 3 of 
10 Galactic low-mass X-ray binaries show absorption of OVII and/or 
OVIII. 



Average n^ hot gas densi ty estimated in the direction of 7 low-mass 
X-ray binaries lYao & Wan j ilOOS,) ). assum ing oxygen abundance rel - 
ative to hydrogen n(0)/n(H) = 4 x IQ-'^ dAnders & Grevesse|[T98^ 
and an ionic fraction of OVII ~ 1 (see, Fig.[3j. 
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= 10' K. r,,= 10' cm 




Fig. 4. Simulated spectra of a photo-ionized plasma associ- 
ated with a NS (blackbody spectrum T=10^ K and lumi- 
nosity of lO^^^ergs"') for different initial hydrogen volume 
densities, located at the distance of r,„„(,r - lO'^^cm. The 
filled dots indic ate the observed optica l /UV fl uxes of RX 
J0720 4-3 1 25 (,Kulkarni & van Kerkwiifl i 1 9981: iKaolan et all 
120031: iMotch et al.ll2003l: lEisenbeiss et al. I l2009h . Fore some pa- 
rameter sets the simulated spectra can explain the power-law 
component with an "optical excess". 



origin of the absorption feature mainly in the ambient medium 
of RX J0720-3125 more likely. 

4. Conclusions 

1. An absorption feature at 0.57 keV is clearly detected 
in the co-added XMM-Newton RGS spectrum of RX 
J0720.4-3125. The detected absorption feature likely is a 
blend. 

2. Most probably, it originates mainly in the ambient medium 
of RX J0720.4-3 125 and may be identified with highly ion- 
ized oxygen (O VII), consisting of Doppler shifted compo- 
nents. Neither an interstellar nor an atmospheric, partial ori- 
gin can be excluded completely. 

3. The observed optical/UV flux excess of RX J0720.4-3125 
compared to the extrapolated X-ray blackbody radiation par- 
tially may be caused by emergent emission of a nearby, rela- 
tively dense photo-ionized cloud. 

Further investigation, in particular simultaneous X-ray and 
optical short-term variations could shed some light on the 
nature of this enigmatic object. 
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spectrum is consistent with the slope of a Rayleigh - Jeans tail 
jvan Kerkwiik & Kulkarnll200ll:lBurwitz et al.ll2003l) . while for 
RX J0720.4-3125 there is a deviation and an additional power - 
law component was necessary to fit the data dKaplan et al.l2003l) . 

In order to estimate the contribution from a photo-ionized 
plasma near an INS, as mentioned above, we have used the 
widespread photoionization code CLOUDY (calculations were 
performed with version 07.02.02 of Cloudy, as described by 
iFerland et al.lll998 ). Despite the restrictions (luminosity of the 
ionizing source, ionization parameter, ie. a combination of loca- 
tion of the cloud relative to the source and spatial density of the 
medium) there are still a number of free parameters for input to 
the CLOUDY code: Covering facto]0 (open or closed geometry), 
chemical abundances, filling factor, turbulence, etc. 

Nevertheless, as shown by our simulations, in principle, the 
existence of a small, relatively dense ambient medium (e.g. riu- 
3.5 X lO^cm"-', covering factor = 0.013, located at the distance 
of 10"^cm from INS, see Fig|4|, may at least partially explain 
the enigmatic "opt i cal excess" observe d at RX J0720.4-3125 
dMotch et all 120031: iKaplan et alj l2003h and provide sufiicient 
OVII ions. 

The satisfactory fit to the observed "optical excess" was not 
unique and may also be provided with other sets of parameters. 

A partial origin of the absorption feature in the local, 
interstellar hot medium (e.g. Local Bubble) cannot be ex- 
cluded completely. However, the non-detection of the absorp- 
tion feature in othe r cases 3 beyond th e Local Bub b le (fo r 
other objects, see, l^ &Wand ( l2005h . iFang et alj (l2006h . 
iBregman & Llovd-DaviesI ( l2007h ) makes the proposition of the 



^ The covering factor is the fraction of 'Xnsr cover ed by gas, as 
viewe d from an ionizing source, ie. £11 4n. For details, see lFerland et alj 
( fl998l) , CLOUDY documents 

* The absorption feature is absen t in the co-added RGS sp ectrum of 
RX J1856.4-3754 (d = 161!]5 pc ( van Kerkwii k & Kaplanl l2007,) . but 
clearly seen in the co-added spectra of RBS1223 and RX J1605.3+3249 
with lower significance (see, http://xmm.esac.esa.int/BiRDJ 
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